t = readmatrix('A.out', 'FileType', 'text');
x = -1:0.001:1;
figure;
for i = 1:6
    yi = t(i, :);
    plot(x, yi);
    hold on;
end
l = {'true cruve', 'N=6', 'N=11', 'N=21', 'N=41', 'N=81'};
legend(l);
saveas(gcf, '../pic/A.eps');

t = readmatrix('A2.out', 'FileType', 'text');
X = t(1,:);
Y = t(2,:);
Y = log(Y);
Y1 = log(1./X./X);
Y2 = log(1./X./X./X);
figure;
plot(X, Y);
hold on;
plot(X, Y1);
plot(X, Y2);
l = {'log max error', 'log N^-2', 'log N^-3'};
legend(l);
saveas(gcf, '../pic/A2.eps');

t = readmatrix('A3.out', 'FileType', 'text');
X = t(1,:);
Y = t(2,:);
figure;
plot(X, Y);
l = {'log max error ratio'};
legend(l);
saveas(gcf, '../pic/A3.eps');
